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ABSTRACT 

Fullerenes have been recently detected in various circumstellar and interstellar environments, raising the question of their formation 
pathway. It has been proposed that they can form at the low densities found in the interstellar medium by the photo-chemical pro¬ 
cessing of large polycyclic aromatic hydrocarbons (PAHs). Following our previous work on the evolution of PAHs in the NGC 7023 
reflection nebula, we evaluate, using photochemical modeling, the possibility that the PAH C66H20 (i.e. circumovalene) can lead to 
the formation of Cso upon irradiation by ultraviolet photons. The chemical pathway involves full dehydrogenation of C66H20, folding 
into a floppy closed cage and shrinking of the cage by loss of C2 units until it reaches the symmetric Ceo molecule. At 10” from the 
illuminating star and with realistic molecular parameters, the model predicts that 100% of C66H20 is converted into Cgo in ~ 10^ years, 
a timescale comparable to the age of the nebula. Shrinking appears to be the kinetically limiting step of the whole process. Hence, 
PAHs larger than C66H20 are unlikely to contribute significantly to the formation of Cgo, while PAHs containing between 60 and 66 C 
atoms should contribute to the formation of Cgo with shorter timescales, and PAHs containing less than 60 C atoms will be destroyed. 
Assuming a classical size distribution for the PAH precursors, our model predicts absolute abundances of C60 are up to several 10“'* 
of the elemental carbon, i.e. less than a percent of the typical interstellar PAH abundance, which is consistent with observational 
studies. According to our model, once formed, Ceo can survive much longer (> 10’ years for radiation fields below Go = 10*) than 
other fullerenes because of the remarkable stability of the Coo molecule at high internal energies. Hence, a natural consequence is that 
Coo is more abundant than other fullerenes in highly irradiated environments. 
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1. Introduction 

The mid-infrared (mid-IR) spectrum of galactic and extragalac- 
tic objects exhibits emission in bands (strongest at 3.3, 6.2, 7.7, 
8.6, and 11.2 pm) attributed to carbonaceous macromolecules, 
more specific ally to polycyclic aromatic hydrocarbons (PAHs, 
lTielensll2008h . In addition to PAH bands, IR emission bands at 
7.0, 8.5,17.4, and 19.0 pm have been reported (ICami et al.l2010l; 
ISellgren et al.ll20foh . and found to match quite cl osely the IR 
active bands of neutral buckminsterfullerene (Ceo, iKroto et ^ 
cm, a cage-like carbon molecule. Additional bands at 6.4, 
7.1, 8.2, and 10.5 pm were observed in the NGC 7023 reflec - 
tion nebula and attributed to the cation (iBerne et al.ll2013l) . 
Carbonaceous macromolecules, especially PAHs, are believed to 
play a fundamental role in the physics and chemistry of the in¬ 
terstellar medium (ISM), and their infrared signatures are com¬ 
monly used as a tracer of physical conditions (especially the 
UV radiation field). PAHs are b elieved to lock up bet ween 5 
and 20% of the elemental carbon (iJoblin & Tielensir201 ll) . while 
Ceo and Cggare found in small abundances in the ISM (at most 
5. 6 X 10 * of the elementa l carbon for the neutral form according 
to lCaste llanos et an20l4 and at most 10“* of the elemental car¬ 
bon for the cation according to lBerne et al.ll201^ . Nevertheless, 
Ceo and CJq are the only species belonging to the family of car¬ 
bonaceous macromolecules which have been specifically iden¬ 
tified in the ISM, and therefore these molecules have attracted 
considerable interest because they open a new possibility to 
probe carbon chemistry and physical conditions in the ISM. 


One question related to fullerenes, and in particu¬ 
lar Ceo, concerns their formation pathway. Recently, “top- 
down” schemes where larger carbon clusters shrink to reach 
Cfio have been proposed dChuvilin et al.l2010l:IZhang et al.l2013l: 

iPietrucci & Andreonill2014l) . and can be opposed to the tradi¬ 

tional “botto m-up” approach where Cm is built up from small 
compounds (iKroto & McKa^ 119881; iHeathI ll992tlHunter et al.l 
1 19941 iDunk et aiT 2013 h. Using inf r ^ed o bservations of the 
NGC 7023 nebula. iBerne & TielensI (12012l) found evidence of 
an increase of the abundance of Ceo with increasing UV field, 
while the abundance of PAHs decreases. This was interpreted by 
these authors as evidence for the formation of Ceo from large 
PAHs (Nc > 60) under UV irradiati on, a top down mech ¬ 
anism similar to the one observed by |Qumlin_et^lJ(|20]^. 
iGarcIa-Hernandez et akl (1201 Ol 1201 ll) and iMicelotta et al.l(l2012l) 
proposed a similar mechanism where the starting materials are 
more complex, such as hydrogenated amorphous carbon instead 
of PAHs. Top-down scenarios are particularly appealing, given 
that the densities prevailing in the ISM are many orders of mag¬ 
nitude too low to allow for a “bottom-up” formation (i.e. start¬ 
ing from small compounds) over reasonable timescale^l In this 
paper, we use a detailed photochemical model coupled to a de- 


* W hile writing the present paper, it was proposed by iPatra et ^ 
1 I2OI4I) that nucleation of C atoms leading to the formation of cages 
can operate in the ISM. However, the densities and relative abundances 
of C and H adopted in this latter work are far from those considered 
here and actually observed in the interstellar medium. 
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scription of the physical conditions in the NGC 7023 nebula. We 
demonstrate that the formation of Ceo by UV processing of large 
PAH molecules is a plausible mechanism to account for astro¬ 
nomical observations. 


cages, but with a sc aling proportional to the number of C atoms 
(I.Toblin et al.lll992t) . The impact of the chosen UV absorption 
cross-sections on the results is discussed in Sect. 14.31 


2. Model for the formation of Ceofrom PAHs 

2. 1. Description of the proposed scenario 


The prop osed scenario of formati on is inspired by the one de¬ 
scribed in iBerne & TielensI (l2012l) and is represented schemati¬ 
cally in Fig.[T] P AHs are assumed to be formed in the envelopes 
of evolved stars ([ Frenklac h & Feigelsonlll98^ICherchneff et al.l 
Il992t iMerino et ^ 2014I) . and then injected in the ISM. Under 
UV irradiation, large PAHs, (Nc > 60) are first fully dehy¬ 
drogenated into small graphene flakes, dehyd rogenation being 
by fa r the dominant dissociation channel (see iMontillaud et al.l 
I2013L and references therein). Further UV irradiation enables 
these flakes to fold into closed cages. Once the cages are 
close d, they can loo se C 2 units if they continue absorbing en¬ 
ergy (llrle et al.ll2006h . Because of the low densities preyailing in 
the considered regions {iih < 10"^cm“^), the reyerse reaction, i.e. 
addition of C 2 is too slow to balance photodissociation. Once 
a system has reached Cgo it will remain in this form for a yery 
long period of time, because of its remarkable stability. 


2.2. The specific case of the transformation of 066^20 into 

Ceo 

In this paper we consider the eyolution of circumoyalene, 
C 66 FI 20 (first molecule in Fig. [T]i, for the physical conditions 
of NGC 7023 (see Sect. [3]l. We restrict ourselyes to a single 
molecule as a starting point, for simplicity, and because mod¬ 
eling a complete population would introduce many free param¬ 
eters that would be irreleyant in this proof of concept study. 
C 66 H 20 is selected as a test case because it gathers the typical 
properties expected for an interstellar PAH (Nc ^ 50, com¬ 
pactness), it can be considered as a potential precursor of Ceo 
and its spectroscopic properties are ayail able in the Cagliari 
Theoretical spectral database of PAH^ (iMalloci et al.l l2007l 
hereafter Cagliari database). The photophysical modeling re¬ 
lies on seyeral key processes: 1) UV absorption which brings 
the molecules to high temperatures, 2) radiatiye cooling, 3) 
unimolecular dissociation, 4) isomerization and more specifi¬ 
cally folding of the graphene flakes and 5) reactions with elec¬ 
trons, H atoms or C^. In the following sections, we present 
the adopted methods to describe these processes which were 
then implemented in t he photochemical model deyeloped by 
IMontillaud et al.l (l2013l) . 


2.2.2. Ionization 

In the frame of our interstellar model, two ionization processes 
should be considered: direct photoionization and thermionic 
emission. Direct ionization occurs after the absorption of a sin¬ 
gle UV photon with an energy aboye the ionization threshold. 
In the case of thermoionic emission, the successiye absorption 
of seyeral UV photons brings the molecule to high temperatures 
and can lead to delayed electron emission. 

For the hydrogenated and dehydrogenated PAHs, the direct 
photoi onization process is included, following IMontillaud et al.l 
(l2013h . For the cages, including Ceo, we do not include direct 
photoionization. This assumption relies on the fact that obser- 
yations indicate that this species is mostly in the neutral form 
in NGC 7023, C^ nbeing localized only yery close to the star 
(iBerne et al.ll2013T) . It cannot be excluded that this is due to an 
efficient recombination of Ceo with electrons from the ambient 
medium, and hence ionization could still be efficient as a sink 
of energy for relaxation. Since the ionization potential of Ceo is 
7.5 eV, only photons with energies 7.5 < hv < 13.6 eV can con¬ 
tribute to ionization. G Nen the spectral profile of the radiation 
held in NGC 7023 (see IMontillaud et'^l20 1 3l) and the fact that 
in this range the mnization yield i s low compared to the total 
photo-absorption (lBerkowitzlll999h . only a small fraction of the 
energy from absorbed UV photons will be lost in photoelectrons, 
and most of it will be conyerted into yibrational energy that can 
eyentually be used for dis sociation. The oth er cages haye sim¬ 
ilar ionization potentials (ISeifert et al.l 1 199^ and are therefore 
assumed to behaye like Ceo ■ Hence, we haye neglected direct 
photoionizaiton for the cages in the photo-chemical model. 

Thermoionic emission b ecomes efficient onl y at internal en¬ 
ergies higher than ~ 30 eV (iHansen et al.ll200ll . It could there¬ 
fore be an efficient relaxation process, competing with photodis¬ 
sociation (Sect. I 2 . 2 . 4 I 1 and cooling by fluorescence from ther¬ 
mally excited electronic states (Sect. 12.2.3l l. As can be seen in 
Fig. |2] the termoionic emission rate becomes larger than fluo¬ 
rescence from thermally excited electronic states only for in¬ 
ternal energies aboye ~ 45 eV. This is well aboye the internal 
energies of dissociation (see Sect. 14.1b and therefore relaxation 
through thermionic emission can be neglected since it will al¬ 
ways be dominated by fluorescence from thermally excited elec¬ 
tronic states. 


2.2.3. Radiative cooling 


2.2.1. UV absorption 

UV absorption is the first step in the photophysical eyolution of 
the considered species. The efficiency of this process in a giyen 
radiation field depends on the UV absorption cross-section. The 
UV absorption cross-section of C 66 H 20 used here is taken from 
the Cagliari database. The same yalue is used for dehydro¬ 
genated PAHs CgeH,, as well as for planar Cee, considering 
that the UV absorption cross-section in a PAH is, in first ap¬ 
proximation, proportional to the number of carbon atoms. For 
Cf,n, w e adopt the UV absorption cross-section from iBerkowitzl 
(II 999 I) . The same absorption cross-section is used for the other 

^ astrochemistry.ca.astro.it/database/ 


Radiatiye cooling can occur through emission of photons in the 
infrared, but when the molecules reach high enough tempera¬ 
tures, it has been shown that fluorescence from thermally excited 
electronic states b ecomes the dominant cooling mechanism (see 
iLeger et aHl 19881 who called this mechanism Poincare fluores¬ 
cence). This process has been obseryed f or internal energies of ~ 
7 eV for the anthracene cation, Ci 4 H|q, (iMartin e t al.ll201 3l) and 
intern al energies of ~ 13 eV for the Ceo anion ( Andersen et al.l 

I2nnih . 


The infrared coolin g rates are c a lculate d using the micro- 
canonical formalism of Lloblin et al.l (l2002ll . This requires the 
full knowledge of the yibrational properties of the considered 
species. For C 66 H 20 , we us e the spectroscopic prop erties from 
the Cagliari database as in IMontillaud et al.l (l2013h . For Ceo , 
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1) Dehydrogenation 



2) Folding 


Fig. 1. Schematic representation of the evolutionary scenario for the formation of fullerenes from PAHs under UV irradiation. 


they are taken from ISchettino et aP (1200 Ih . For dehydrogenated 
PAHs CgeH,, and cages, we use the DFT vibrational frequencies 
of C 66 H 20 from the Cagliari database, after removing the C-H 
vibrational modes and 3 modes per missing C atom. The cooling 
rates calculated with this approach for Ceo are reported in Fig. |2] 
Cooling by fluorescence from thermally excited elec- 
tronic states is calcul ated using the formalism described in 
IChunka & KlotsI (Il997h . The cooling rate is given by: 

k{T) = (2clao)Y^fi{hvilmc^f expi-hvilksT) (1) 

where ao is the Bohr radius, fi and v, are respectively the oscilla¬ 
tor strengths and frequencies of the electronic transitions and T 
is the temperature of the molecule. For C 66 H 20 and for the dehy¬ 
drogenated PAHs we use the v,- and fi values for C 66 H 20 from the 
Cagliari database (Table [T]i. For the cag es, we use the energies 
and oscillator strengths of Ceo taken from IChupka & Klotsll 19971 
reproduced in Table [1] Note that these oscillator strengths are a 
convenient approximation to calculate the radi ative cooling but 
are no t true oscillator strengths. Nevertheless, IChupka & KlotsI 
(Il997h have shown that this approximation is in very good agree¬ 
ment with rate calculations including detailed molecular prop¬ 
erty information. Using the microcanonical formalism described 
above for the calculation of the IR cooling rates, we derived the 
relation between internal energy E and temperature T for the dif¬ 
ferent species. Using Eq.[Tl this allowed us to derive the radiative 
cooling rates as a function of E. For the case of Ceo these rates 
are shown in Fi g. |2]l and are i n goo d agreement with the earlier 
work models of lT^ita et al.l (l2003h . 

2.2.4. Photocdissociation: dehycdrogenation and shrinking by 
ioss of C 2 

Photodissociation is treated using the statistical approach, and 
based on the inverse Laplace transform of the Arrhenius 
equation (see details in iMontillaud et al.l 1201 3l and references 
therein). The parameters involved in this calculation are the ac¬ 
tivation energy Eq, the pre-exponential factor A, and the vi¬ 
brational densities of states that are computed using the vi¬ 
brational frequencies. For PAH dehy drogenation (first step i n 
Fig . [U, we use the parameters given in IMontillaud et all (l2013h . 
For the shrinking of cages by loss of C 2 molecules (step 3 in 
Fig. [D, experimental and theoretical results are scarce, except 
in the case of Cgo and Ceo ^ for w hich published mo lecular data 
have a large scatter, as noticed by iMatt et akl (1200 ll) . These au¬ 
thors reanalyzed published experiments of Ceo ^ dissociation and 


Table 1. Frequencies and oscillator strengths adopted in the ap¬ 
proximation to calculate the cooling by fluorescence from ther¬ 
mally excited electronic states (see text for details). 


A (A) y, (Hz) f, 

_ C66H20 _ 

3470 5.48 X 10“* 0 I 2 

4980 6.02 X 10'^ 0.24 


Cso 


4540 

6.60 X 10" 

0.214 

3485 

8.60 X 10" 

1.17 

2725 

1.10 X 10" 

2.02 

2240 

1.34 X 10" 

1.30 


showed that when using similar and consistent sets of molecu¬ 
lar data, most experiments lead to similar results, with a pre¬ 
exponential factor of 5 X lO^^s ' and a dissociation energy of 
10.0 + 0.2 eV. The activation energy should be equal to the 
dissociation energy if there is no barrier for the reve rse reac¬ 
tion, i.e. C 2 addition. Experime ntal (iLaskin et al.lll99'^ and the¬ 
oretical (lYi & Bernholclll993l) studies indeed indicate the ab¬ 
sence of such a bar rier. More recent studies are generally in line 
with these results (iDiaz-Tendero et aT]l2003^ iGluch et'Sl2004 
among others). We consider each shrinking step from Cee to C 58 
individually (Eqs. 2 to 6 in Ta bleO. Dissociation energies have 
been determined theoretically dZhang et al.lll992li and the value s 
differ significantly from experimental data (iGluch et al.ll2004l) . 
The theoretical values were computed directly for neutral cages 
in a consistent manner for all species. In contrast, the experimen¬ 
tal values were derived from measurements on fullerene cations. 
We used the theoretical values in our standard model (Table |2]i. 
This choice has the advantage of providing consistent energies 
for all cages in this study, but with a r elatively high val ue for 
Ceo (11.1 eV), compared to the results of iMatt et aklfcoOlli (10.0 
eV). To reconcile our ch oice of binding energies with the re¬ 
sults of iMatt et al.l (1200 Ih . we adopted a pre-exponential factor 
of 2 X ipZQs-i which lea ds to a dissociation rate comparable to 
that of iMatt et~^ (1 19991) (their Fig. 8 ). We used the same factor 
for the C 2 loss of all cages. The impact of all these assumptions 
is discussed in Sect. 14.31 
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Fig. 2. Cooling rates for Ceo as a function of the internal en¬ 
ergy of the molecule. The infrared emission and cooling by fluo¬ 
rescence from thermally excited electronic states are calculated 
(labelled Poincare in the figure) following the formalism de- 
scrib ed in this paper. Th e thermionic emission rates are taken 
from iHansen et alJ (l2003h . 


2.2.5. Folding 

For the folding of the carbon flakes (step 2 in Fig.[T]), there are, to 
our knowledge, no experimental data to constrai n the activation 
energy and pre-exponential factors. Recently, iLebedeva et alJ 
( 12012 h conducted molecular dynamics simulations of the folding 
of the C 96 and C 384 graphene nanoflakes. For these two systems, 
they ran simulations at several temperatures and found that the 
rates follow an Arrhenius law from which they were able to de¬ 
rive the effective activation energies and pre-exponential factors. 
They show that these parameters do not vary significantly with 
sizqj hence we adopt their values reported for C 96 to describe 
the folding (see Table|2]i. 

2.2.6. Reaction with electrons, H atoms or 

The other processes of interest include the recombination with 
electrons and the reaction with H atoms. During the dehy- 
drogenation step all th ese processes are included, following 
iMontillaud et alJ (l2013h . We do not include the addition of C 2 , 
because the abundance of this species is expected to be very low 
in the cavity of NGC 7023. The cation is however abun¬ 
dant, and could react efficiently with cages. There are to our 
knowledge no studies of the kinetics of this reaction. We es¬ 
timated the maximum contribution of this process in the case 
of C fin using the Lang evin rate. A polarizability of ~ 80 A^ for 
Ceo ^ZoDe et al.ll2008ll leads to ^Langevin ~ 6 x 10“® cm^s ', and 
assuming a high abundance of of 3 x 10 «h, one gets an ef¬ 

fective addition rate of ~ 2 x 10“^^ «h (that is < 4 x 10“^ 
s“' for hh < 2x lO^cm^^). This rate must be compared to the rate 
of the reverse process, i.e. photodissociation of Cgj. We mod¬ 
eled this latter process with the same formalis m as described 
in Sect. [3] using a dissociation energy of 4 eV (ISlanina & Lee! 

^ It should be noted that this assertion concerns the 100-400 C atoms 
range and may not extend to lower C numbers, yet as discussed in 
Sect. 14.31 the activation energy would need to be different from the 
adopted one by a large factor to alter our results 


Table 2. Activation energy and pre-exponential factors for the 
key reactions considered in the model 





Act. energy 
E, (eV) 

Pre-exp. factor 
A(s-‘) 

(1) 

^planar 

- 

4.0 

1 x 10'^ 

(2) 

^cage 

^66 

cz + C 2 

8.1“/9.4*’ 

2 x lO^® 

(3) 

(-^cage 

'^64 

*-62 + *-2 

9.4“ / 9.0* 

2 x 10^° 

(4) 

(-^cage 

'^62 

r‘^age ^ 

*-60 + *-2 

6.0“/8.1* 

2 x lO^® 

(5) 

'^60 

r^cage ^ 

*-58 + *-2 

11.1“/ 11.2* 

2 x 10^® 

(6) 

^cage 

'^58 

f~,cage ^ 

^6 + ^2 

8.7“/9.7* 

2 x 10^® 


Notes, a. Using acti vation energ i es of IZhang et al.l (Il992h : b. Using 
activation energies of iGluch et all ( l2004h . 


le-2 2e-2 3e-2 4e-2 5e-2 



Fig. 3 . Physical conditions in NGC 7023. Density profile de¬ 
rived from observations of the far-infrared dust emission in the 
cavity of NGC 7023 (blue with error-bars) and exponential fit 
to this profile used in the photochemical model (red line). The 
black line shows the adopted profile for the intensity of the radi¬ 
ation held Gq. 


1994 ) and a pre-exponential factor of 1.6 x 10^^ s“^ (iKlotsI 
199 ih . We And photodissociation rates from ~ 7 x 10“^ s”' to 
~ 2x10“^ s ' for the various astrophysical conditions considered 
in Sect. 14.21 Therefore, the photodissociation of Cei is estimated 
to be generally much faster than C^ addition in the regions of 
interest, and we do not include this latter process in our study. 

3. Astrophysical environment: NGC 7023 


Since NGC 7023 is the t emplate interstellar source for the study 
of fu l lerene formation ( Sellgren et alJ 1201 Ot iBerne & TielensI 
1201 21: IMontillaud et al JI2013l) . we will model the photochemi¬ 
cal evolution of C 66 H 20 for conditions found in this nebula. The 
formation of Ceo is thought to occur inside the cavity between 
the star FID 200775 and the PDR situated 40” at the North- 
West of the star. Here, we study five positions at 5, 10, 15, 20 
and 25” from the star , situated on a cross-cut (see Fig. 1 in 
IMontillaud et alJl2013h that goes from the star to the PDR. On 
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lO'^ 10'^ 10'^ 10'^ 10° 10^ 10^ 10^ 10'* 10^ 10° 10^ 


time [ yr ] 


20 , 

Fig. 5. Model results for the time-evolution of the abundance of PAHs ( 2 CeeHn), and cages (Cgg Cgg**', Cg 4 *‘',Cg 2 ^‘', 

n=\ 

and C 2 . These are normalized to the abundance of C 66 H 20 at f = 0. The model calculations are done for distances of 5, 10, 15, 20 and 
25” from the star HD 200775 in NGC 7023. Note that the dissociation of C 2 is not treated in our model, hence the C 2 abundances 
reported are only provided as an indication of how much C 2 is formed from the dissociation of cages. Similarly, the dissociation of 
the is not included in the model and therefore its abundance is provided as an indication of how much Ceo is being destroyed. 
The gray shadowed regions indicate the uncertainty range implied by uncertainties of a factor of 10 on the rates of C 2 loss (see 
Sect. 14.31 for details). The approximate age of NGC 7023 is given, with orange shadowed regions representing the uncertainty on 
this value (see Sect.[3]for details). 


the same cross-cut, we derive the spectral energy distribution of 
dust emissi on using archival da ta from the Herschel space tele¬ 
scope Isee lAbergel et alJl2010l) . Using these SEDs, we derive 
the column density on the line of sight using a value for dust 
opacity at 250 u ra. of 1.14xl0~^^cm~ per H atom (see details 
on the method in iPlanck Collaboration XX\1l201 ih . This yields 
a column density pr ofile, which we div ide by the thickness of 
the region derived bv iPilleri et alJ (1201 2h based on the modeling 
of the PAH emission profile, i.e. 0.13 pc. The resulting density 
profile is shown in Fig. [3 To obtain an analytical description of 
this profile, we fit it with an exponential law (see Fig. O. The 
intensities of the radiati on field at 5, 10, 15, 20 and 25” from 
the star are derived as in iMontillaud et all (1201 3h and the corre¬ 
sponding profile of Go (UV field intensity in terms of the Habing 
field, which corresponds to 1.6 x 10 ~^ erg cm~ ^ s ' when inte¬ 
grated between 91.2 and 240 nm. lHabin^ 1968h as a function of 


distance to the star is presented in Fig. [3 In order to derive the 
temperature of the gas at these positions PDR modeling is re¬ 
quired. For the conditions described above the gas temperature 
derived by the Meudon PDR model typically ranges between a 
few 100 and 1000 K. Here we adopt a characteristic tempera¬ 
ture of 300 K for the photochemical model but the exact value is 
not critical (see Sect. l4.31 l. Finally, a parameter that is key when 
comparing the results of the model with observations is the age 
of the nebula, i.e. the time from when the cloud started to receive 
UV photons. This number cannot be derived directly and hence 
other indicators have to be used. One is the age of the illumi¬ 
nating star HD 200775. It is difficult to give a precise “age” for 
such a young and massive star, ho wever, it is most likely rang ing 


between lO'' and 5 x 10^ years (AIecian_e^^ 200^ 12013 and 
Alecian 2014 priv. comm.!. iBerne & Tielens ( 2012h adopted a 
value for the age of the nebula ranging from a few lO"' to 1 x 10 ^ 
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E(eV) 

Fig. 4. Probability density function of dissociation of the 
Ceo molecule as a function of internal energy, at a distance of 
5” from the star (Go ~ 2 x 10^). Energies Epeak at maximum 
dissociation rate {D(E) = kdissiE) x n{E)), Eio„ and Eup (where 
D{Eup) = D{Eio„) = D(Epeak)I^OQ and E^p > Eio„) are de¬ 
picted. The values of Epeak, Eio„ and Eup for all cages are given 
in Table [3] 


years. The analysis of the dynamical properties of the warm neu¬ 
tral gas tr aced by the lin e with Herschel indicate an age of ~ 
0.5 Myrs (iBerne et alJl2012h . However, more recent calculations 
focusing on the photo-evaporation flow in the NGC 7023 North 
PDR indicate ages as low as 1.6 x 10"^ years (Berne et al. in prep). 
We will use 10^ years as the reference age, but we will consider 
a range between 10^ and 5 x 10^ as an age uncertainty range. 

4. Results 

4.1. Internal energies 

In order to understand the processes at play in the photochem¬ 
istry, it is useful to evaluate the probability density functions 
(PDFs) of the internal energies of the species that photodisso- 
ciate with the model. Fig. |4] shows such a PDF for Ceo at a dis¬ 
tance of 5” from the star. In Table [3] we report the main charac¬ 
teristics of these PDFs for all the cages considered in the model: 
the peak energy Epeak i-S- the energy of maximum dissociation 
rate D{E) expressed in s ', and the lower and upper bounds Eiow 
and Eup defined by D{Eup) = D{Eia„) = D{Epeak)!lOQ with 
Eup > Eiow (see Fig.|4]for a graphical definition of these param- 
et ers). These values a re calculated using the activatio n energies 
of IZhang et al.l(ll992ll . and those of iGluch et al.l(l2004ll (Table|2]i. 
The results presented in Table |3] indicate that the internal ener¬ 
gies required to dissociate the molecules are high, i.e. > 15 eV. 
Photons in the NGC 7023 nebula have a maximum energy set 
by the Fyman limit, i.e. 13.6 eV. Hence, the molecules that dis¬ 
sociate must have absorbed multiple photons, and therefore the 
photochemistry of these species is completely governed by these 
multiple-photon absorptions. In the specific case of Ceo, the en¬ 
ergy needed is at least 27 eV, i.e. requiring the absorption of at 
least three photons in order to dissociate. 


Table 3. Properties of the probability density functions of dis¬ 
sociation of cages (see example of such a function in Fig. Ell : 
Epeak, Eiow and Eup. All values in eV and derived for a distance 
of 5” from the star (Go ~ 2 x 10^). 




a. 



b. 


Species* 

^low 

^peak 

Eup 

^iow 

^peak 

Eup 

C66 

22 

32 

39 

26 

35 

43 

C64 

24 

35 

43 

24 

34 

41 

C62 

15 

20 

28 

22 

30 

36 

Cso 

27 

36 

47 

27 

36 

47 

Css 

22 

30 

37 

24 

33 

41 


Notes. *Cages only are pre sented in this table, a. Using act ivation ener- 

! :ies ofl^ang et al.llll992li : b. Using activation energies of iGluch et all 
200 ^ . 


4.2. Abundances 

The results of the model are given in Fig. |5] which presents the 

20 

time-evolution of the abundance of PAHs (2 CeeH,,), and cages 

(CPianar^ C“/^ C“"^C“*^ C“*^ C™/'') andQ, at distances of 5, 
10, 15, 20 and 25” from the star HD 200775 in NGC 7023. For 
these five positions, full dehydrogenation occurs very quickly (a 
few 10s of years at 25” and a few days at 5” from the star). This 
implies that C 66 H 20 is quickly destroyed in NGC 7023, however 
larger PAHs could survive over longer timescales. Once the pla¬ 
nar Cee is formed, it immediately folds and forms a cage. This 
implies that graphene flakes are only transient species and there¬ 
fore unlikely to be detected in the ISM. The Cm cage, instead, 
can survive for a relatively long period of time (few 10s of years 
at 5” from the star and up to several 10^ years at 25” from the 
star). It is therefore the shrinking step, where cages loose C 2 
units, that limits the efficiency of the Ceo formation process in 
this model. Since each of the C 2 loss process is time consum¬ 
ing, it appears unlikely that PAHs with sizes Nq > 66 C atoms 
will contribute significantly to the formation of Ceo ■ At distances 
shorter than 10” from the star the cages can shrink to Ceo and af¬ 
ter 10^ years all the C 6 eH 2 o has been converted to Ceo ■ Once the 
molecules have reached Ceo, it takes a very long time to destroy 
them: Ceo survives for lO’s of Myrs at distance larger than 10” 
to the star (radiation fields below Go of a few 10^). Only at the 
closest position to the star (5”, Go ~ 2 x 10^) is Ceo destroyed 
efficiently and hence it is predicted that its abundance will de¬ 
crease after a few 10"^ years. 

4.3. Sensitivity of model results 

Is the model very sensitive to the adopted parameters, and if so, 
which are the critical parameters ? 

First, we have checked that the dehydrogenation timescale 
is always much smaller than the cage formation plus shrinking 
timescale for the physical conditions of the considered positions. 
Therefore, it is equivalent to start our calculations with the pla¬ 
nar Cee rather than with C 66 H 20 . All the following tests have 
been conducted with the above assumption, except for the gas 
temperature test. We first test the choice in the vibrational fre¬ 
quencies for the system. We And the exact values of the fre¬ 
quencies to have a negligible impact on the results, as long as 
the number of vibrational degrees of freedom and the orders of 
magnitude of the frequencies are correct. We have varied the 
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gas temperature in the model from 100, to 300 and 1000 K 
and found that this has a negligible impact on the Ceo formation 
timescale. The UV-visible absorption cross section of the cages, 
which is not well known, has a somewhat larger impact on the 
results. For instance, if for the cages we u se the cross section 
per C atom of CfifiH^n (Malloci et al.ll2007h instead of the one 
of Ceo (iBerkowitzl 1 1 9991) 7 we obtain a variation by a factor of 
three in the timescale of formation of Ceo- Multiplying the pre¬ 
exponential factors of C 2 loss by 0.1 or 10 simultaneously for 
all cages we observe an inversely proportional variation of the 
Ceo formation timescale. This effect in the abundance of Ceo is 
illustrated by the shadowed region in Fig.|5] When varying indi¬ 
vidually these coefficients, we hnd that the results are impacted 
mainly by the value for Cge —» C 64 - 1 - C 2 , and marginally for 
C 64 ^ C 62 + C 2 , while changes for the two other shrinking reac¬ 
tions do not impact the results. We comp ared our standard r esults 
(i.e. using the activation energies from IZhang et al.lll992h with 
the results obtained when usi ng the experimental dissociation 
energies of iGluch et al.l (l2004ll . The Ceo formation timescales is 
slightly increased, mainly due to the higher stability of 
in the experimental data set. Overall, the effect does not affect 
the comparison with observations and discussion which follow. 
Flence from now on we only cons ider the results of the model 
using the activation energies from IZhang et alj[l992i The main 
source of uncertainty appears to arise from the pre-exponetial 
factors chosen for the shrinking steps. We will therefore con¬ 
sider their effect in the comparison with observations. 


5. Comparison with observations 


In this section we compare the results of our model with the 
observation of Ceo formation in NGC 7023. To be completely 
accurate, this comparison would require that the model contains 
a whole size distribution of PAFIs in agreement with observa¬ 
tions. This requirement is beyond the scope of this paper, but the 
agreement between the observations and the model can be tested 
bearing this limitation in mind. 

The comparison is primarily based on confronting the max¬ 
imal abundance of Ceo observed in NGC 7023, at a distance of 
11” from the star, with the value obtained in the model for a dis¬ 
tance of 10”. Taking into account the uncertainty on the age of 
NGC 7023 and the uncertainty on the pre-exponential factors, 
the model-predicted abundances of Ceo ranges between a mini¬ 
mum value of 1.3 X 10“^ and 0.98, relative to C 66 H 20 (Fig.|5]l. To 
convert this into an absolute abundance of Ceo one needs to know 
the abundance of C 6 eH 2 o relative to the total PAH population, 
and to know the t otal abundance PAHs in NGC 7023. The latter 
was estimated bv iBerne & Tielen^ (1201 2l) to be - lx 10“^ 

expressed in fraction of carbon locked in PAHs. The maximum 
abundance of C 6 eH 2 o relati ve to the to t al PAH population, a, 
can be estimated fo Uowing Pillm_gt_^ (l2009l) . using the PAH 
size distribution of iDesert et al.l ( ToOOl) . which yields a value 
a - 5.3 X 10“^. Thus, the fraction of carbon locked in C 66 H 20 
is 5.3 X 10“^ X 7 X 10“^ = 3.7 x 10“"^. With this value, the frac¬ 
tion of carbon locked in Ceo predicted by the model at a dis¬ 
tance of 10” from the star ranges between /5“ = 4.8 x 10“® and 


- 3.6 X 10 (where the lower and upper limit of the range 
include the uncertainty on the pre-exponential factor and the age 
of the nebula). This range of values is in agreement with the one 
derived bv iBerne & Tielen^ (1201 2h of = 1.7 x 10 So far 

we have considered C 66 H 20 as the only precursor of Ceo ■ As 
mentioned in Sec. 14.21 it is reasonable to expect that only PAHs 
bearing between 60 and 66 C-atoms will be able to form Ceo, 


and hence can increase the hnal abundance of Ceo- Hence, it is 
realistic to calculate a incorporating all the species bearing be¬ 
tween 60 and ~66 C atoms, which yields or = 3.8 x 10“^. With 
this value, the abundance of Ceo at 10” from the star is predicted 
to range between /^“ = 3.5 x 10“^ and - 2.6 x 10~^ 
in good agreement with the value derived by iBerne & TielensI 
(12012h of = 1.7 X 10 The comparison can be extended 
to the other positions (15, 20 and 25”) where the abundance of 
Ceo has been measured. At these three positions, the model de¬ 
rived ranges of abundance are also in good agreement with the 
observed abundances. However, as the distance grows the range 
of acceptable values predicted by the model become broader and 
less constraining. Overall, the present comparison demonstrates 
that our scenario is compatible with observations within the un¬ 
certainties on the molecular parameters and the age of the neb¬ 
ula. 


6. Discussion 


6.1. Comparison to other models 

Models o f Cm formation in a top -dow n mechanism w e re pre - 
sented in iBerne & TielensI (l2012h and iMicelotta et alJ (l2012h . 
The scenarios in t hese papers are quite sim ilar, the main dif¬ 
ference being that IBerne & Tielens (|2Q12|) consi d er PA Hs as 
the starting ingredient while for Micelotta et al.l (l2012l) it is 
nanometer-sized “arophatic” clusters. In both of these studies 
the evolution of the species is described by a “thermal model” 
(Arrhenius equation), where the driving parameter is the activa¬ 
tion en ergy. For the los s of C atoms at the edge of graphene 
sheets, IBerne & TielensI (l2012h used a value tuned to 4.5 eV 
to obtain reasonable formation efficiencies. Yet, as noted by 
IMicelotta et alJ (l2012h this value is somewhat arbitrary and the 
obtained Ceo formation efficiencies rema in low. To explain the 
efficient formation of Ceo in evolved stars. IMicelotta et ^(l2012l) 
used the results fr om the molecular dynamics simulations of 
IZheng et"^ (l2007ll . and extracted a unique activation energy of 
0.35 eV for the shrinking reaction by loss of C 2 . However, this 
value is a factor of ~ 20 smaller than what is generally accepted 
for this reaction (see Table 1), leading to over-estimations of the 
rates by many orders of magnitude, ca sting serious doubts on the 
conclusions of IMicelotta et al.l (l2012ll . In summary, since these 
simple models do not allow the molecules to reach high inter¬ 
nal energies through multiple photon absorptions, high forma¬ 
tion rates of Ceo require activation energies tuned to values that 
are too low to be physical. 

In summary, since these simple models do not consider mul¬ 
tiple photon absorptions, the molecules never reach the high 
internal energies necessary to a high formation rate of Ceo ■ 
Therefore, they strongly underestimate the Ceo formation yields, 
unless one assumes activation energies that are too low to be 
physical. Instead, the approach presented here, including a com¬ 
plete description of the photochemical processes, allows, using 
realistic molecular parameters, to predict formation efficiencies 
that are in agreement with observations. 


6.2. Comparison to experimental results of gas-phase cage 
formation 

After we submitted the present paner. l^en et alJ (1201 4l) reported 
experimental results in which they demonstrate that Ceo can 
be formed in the gas phase through laser irradiation of larger 
PAH molecules. These results support the idea put forward here 
that Ceo can be formed following a top-down scheme. However, 
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IZhen et ^ (l2014t) suggest that in some cases the conversion of 
graphene into cages could involve a prior step of C2 loss. This is 
in contradiction with our photochemical model in which the dis¬ 
sociation by loss of C2 is less efficient than the folding by many 
orders of magnitude, and therefore folding always precedes the 
loss of C2. This is mainly because the activation energy for fold¬ 
ing (which relies on molecular parameter s derived from molec- 
ular dynamics simulations performed by iLebedeva et al.ll201^ 
is a factor of ~ 2 smaller than the activation energy for C2 loss. 
Further experimental investigation is necessary, in particular to 
quantify the activation energies involved in the folding and C2 
loss by flakes, so that they can be included in our model. 

6.3. Stability of Ca) 

Table [3] demonstrates that Ceo can reach particularly high in¬ 
ternal energies. In this state, fluorescence from thermally ex¬ 
cited electronic states becomes a very efficient cooling chan¬ 
nel, postponing the dissociation and conferring an increased 
stability to CeoD Hence, in the frame of the top-down mech¬ 
anism detailed here and because of its remarkable stability, 
Ceo is expected to be the most abundant of all fullerenes in 
highly UV irradiated environments in space. It is interesting to 
note t hat similar kinetic sta bility arguments have been put for¬ 
ward (iFedorov et alJ I20T1!) to explain the high abundance of 
Ceo rel ative to other fullerenes observed in labo ratory experi¬ 
ments (iKroto et al.lll985UKratschmer et al.lll990h . As discussed 
in Sect. 14.21 Ceo is expected to survive for 10s of Myrs in intense 
radiation fields. In the diffuse ISM, where the radiation field is 
several orders of magnitude smaller, Ceo is therefore expected to 
resist to radiation for even larger timescales. In such conditions, 
other energetic processes must be invoked to destroy Ceo, such 
as shocks or cosmic rays. 

6.4. Isomerization of Ceo 

There exists a large number of Ceo isomers, and hence there 
is a-priori no reason for the icosahedral //r-Ceoto be the one 
present in the ISM. /h-Ceohas the lowest energy amongst all 
Ceo iso mers and the closest isomers lie about 2 eV higher in 
energy dRaghavachari & Rohlfingiri992h . Upon UV absorption, 
the isomerization reaction from any Ceo isom er towards Ih- 
Ceo re quires only 5.4 eV of activation energy (lYi & Bernholcl 
Il992h . and is therefore expected be fast in the conditions in¬ 
vestigated here. Conversely, the isomerization reaction from Ih- 
Ceo to another Ceo isomer requires ~ 5.4 -H 2 = 7.4 eV, and will 
be much slower than the previous reaction. Therefore, relatively 
quickly, the population of Ceo molecules will be dominated by 
the icosahedral isomer. 

7. Conclusions 

We have presented the first detailed photochemical modeling of 
Ceo formation from PAHs in space following a top down scheme. 
This model is calculated for a single molecule (CeeH2o) as a 
starting point, and the key processes are explicitly described: 
UV photon absorption (including multiple photon events), ra¬ 
diative cooling and dissociation. Because the involved activa- 

^ Note that since multiple photon absorptions are rare, most 
Ceo molecules in the ISM are situated at much lower internal energies 
and the population-averaged emission at a given time remains domi¬ 
nated by infrared transitions. The visible emission of Cgo in the ISM is 
therefore probably very weak. 


tion energies in this top-down mechanism are high (especially 
for the shrinking steps), multiple photon absorptions dominate 
the photochemistry. Using the physical conditions which pre¬ 
vail in the NGC 7023 reflection nebula, we find that Ceo can be 
formed from C66H20 over a timescale of about 10^ years with a 
high efficiency in highly irradiated regions. Assuming that only 
PAHs containing between 60 and 66 C atoms are precursors of 
Ceo and a classical size distribution for PAHs, the comparison 
between the modeled and observed abundances of Ceo shows a 
good agreement, within the uncertainties of the model. These un¬ 
certainties can be reduced once a better characterization of the 
reaction rates of the shrinking of cages is available. While devel¬ 
oped for the physical conditions of NGC 7023 as representative 
of the interstellar medium, the scenario and model described in 
this paper could also be applied to the irradiated circumstellar 
gas of planetary nebulae, if the physical conditions in these en¬ 
vironments can be characterized in detail. 
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